Psychophysiology, (2012), Wiley Periodicals, Inc. Printed in the USA. 
Copyright © 2012 Society for Psychophysiological Research 
DOI: 10. 1 1 1 1/j. 1469-8986.201 1 .01341.x 



Stability of heart rate variability indices reflecting 
parasympathetic activity 



KATJA BERTSCH, ab DIRK HAGEMANN, C EWALD NAUMANN, b HARTMUT SCHACHINGER, d and 
ANDRE SCHULZ d 

"Department of General Psychiatry, Center for Psychosocial Medicine, University Hospital Heidelberg, Heidelberg, Germany 
b Department of Psychology, University of Trier, Trier, Germany 

c Section of Personality Research, Institute of Psychology, University of Heidelberg, Heidelberg, Germany 
d Institute of Psychobiology, University of Trier, Trier, Germany 



Abstract 

Heart rate variability (HRV) is a measure of autonomic influences on heart rate that has frequently been used as a 
transsituationally consistent biomarker for cardiovascular health and emotional or cognitive functions. The psychometric 
properties of HRV however remain unclear. In the present study, we examined the reliability and temporal stability of 
parasympathetic HRV measures and estimated the portion of variance explained by transsituationally consistent trait 
variance and by effects of the situation and person-situation interaction with structural equation modeling. The results 
show good reliability of indices reflecting central parasympathetic control over heart rate and that about 40% of the 
variance of a single HRV measurement can be explained by effects of the situation and person-situation interaction. An 
aggregation across at least two measurements may be recommended when using HRV as a transsituationally consistent 
biomarker or trait. 

Descriptors: Heart rate variability, Electrocardiogram, Latent— state trait model, Autonomous nervous system, Para- 
sympathetic nervous system, Respiratory sinus arrhythmia, RSA 



Heart rate variability (HRV) is a noninvasive, simple, and fre- 
quently used measure of autonomic influences on heart rate. 
Evidence indicates that HRV indices of sympathetic and parasym- 
pathetic activation pattern reflect biomarkers not only for cardio- 
vascular health, but also for complex patterns of brain activations 
(e.g., Aasman, Mulder, & Mulder, 1987; Hansen, lohnsen, & 
Thayer, 2003; Sloan, Korten, & Myers, 1991). HRV has also been 
used as an objective descriptor of individual differences in regulat- 
ing emotional conditions or cognitive functions (e.g., Ahs, Sollers, 
Furmark, Fredrikson, & Thayer, 2009; Hansen et al,, 2003; Kemp 
et al., 2010; Thayer, Hansen, Saus-Rose, & Johnsen, 2009; Thayer, 
Yamamoto, & Brosschot, 2010). HRV thus provides information 
on the nervous system's ability to organize homeostatic responses 
in adaptation to the physiological and psychological requirements 
of the present situation (e.g., Althaus, Mulder, Mulder, Van Roon, 
& Minderaa, 1998; Appelhans & Luecken, 2008; Hjortskov et al., 
2004; Segerstrom & Nes, 2007). In other words, HRV is interpreted 
as a transsituationally consistent disposition or trait that reflects a 
person's adaptability to situations (e.g., Thayer et al., 2009). In 
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most of the studies, HRV is assessed by a single short-term 
electrocardiographic (ECG) measurement in a resting or paced 
breathing condition or during the performance of a cognitive or 
affective task (see, e.g., Sandercock, Bromley, & Brodie, 2005). 

The psychometric properties, that is, the transsituational con- 
sistency (i.e., consistency across differences measurement occa- 
sions and experimental conditions) and temporal stability of 
HRV measures, however, remain unclear (Sandercock, 2007). 
Several researchers have investigated the stability of HRV meas- 
ures across several measurement occasions at rest or during cog- 
nitive or physiological challenges by calculations of test-retest 
correlations (r), intraclass correlations (ICC), or coefficients of 
variation (CV; e.g., Kobayashi, 2009; Lord et al., 2001; Pinna 
et al., 2007; Pitzalis et al., 1996; Schroeder et al., 2004; Sinnreich, 
Kark, Friedlander, Sapoznikow, & Luria, 1998). A meta-analytic 
review (Sandercock et al., 2005) revealed inconsistent results 
regarding the effects of elapsed time between measurement occa- 
sions on stability of HRV measures, ranging from poor stability 
for resting measures with very short test-retest intervals (e.g., 
CV = 45%; Lord et al., 2001) to good stability with test-retest 
intervals of several months (e.g., spontaneous breathing: .80 £ 
ICC < .87; paced breathing: .71 < ICC < .88; Kobayashi, 2009; 
spontaneous breathing: .65 < ICC < .77, 6% < CV < 12.1%, 
paced breathing: .64 < ICC < .82, 6.1% < CV < 10.7%, Sinnre- 
ich et al., 1998; spontaneous breathing: .76 £ r £ .86; cognitive 
challenge: .27^r<\38; Sloan, Shapiro, Bagiella, Gorman, & 
Bigger, 1995). These inconsistent findings are difficult to interpret 
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because a low temporal stability may be due to measurement errors 
or due to a true change across the time interval. This true change 
may, in turn, be due to systematic effects of the measurement 
situation and the interaction between the situation and the person. 
Thus, the finding of a rather low temporal stability of the measure 
may suggest an important role of situation-specific state variability 
in HRV and questions the use of a single HRV measurement as a 
stable biological marker or trait. The current study therefore esti- 
mated the quantitative proportions of situation-specific effects and 
consistent physiological trait factors in HRV measurement. Fur- 
thermore, as several researchers pointed out before the importance 
of multiple HRV assessments, we aimed to explore how many 
measurement occasions are required for conclusions about HRV 
measures as trait variables. We focused on HRV measures reflect- 
ing parasympathetic activity because (a) different parameters in the 
time and frequency domain are available to estimate central para- 
sympathetic control over heart rate and (b) activity in the high- 
frequency band is considered to clearly represent vagal tone, 
whereas the low-frequency band reflects both sympathetic and 
parasympathetic activity and may also be affected by other factors, 
such as diurnal variations (cf. Berntson et al., 1997; Task Force of 
the European Society of Cardiology and the North American 
Society of Pacing and Electrophysiology, 1996). 

Therefore, the purposes of the present study were (1) to assess 
the reliability of HRV measures reflecting parasympathetic activity 
in three different tasks and (2) to estimate the portion of parasym- 
pathetic HRV variance explained by transsituationally consistent 
trait variance and by the situation and person-situation interaction 
with structural equation modeling (SEM). In particular, we used 
latent state-trait (LST) theory (Steyer, Ferring, & Schmitt, 1992; 
Steyer & Schmitt, 1990), an extension of classical test theory, 
which takes into account that no measurement takes place in a 
situational vacuum. Whereas in classical test theory the observed 
variables (e.g., HRV) are decomposed into a true score variable and 
measurement errors, in LST theory the decomposition is performed 
in two steps: In a first step, the observed variables within each 
situation of measurement are decomposed into a latent (not 
observed) state variable and into measurement errors. In a second 
step, all latent state variables are decomposed into a common latent 
trait variable and into occasion-specific state residuals. These 
residuals comprise the effects of the situation and the interaction 
effects of person and situations on the latent state. The two-step 
decomposition of the observed variables thus allows the distinction 
between transsituationally consistent differences (i.e., the latent 
trait), situation-specific fluctuations (i.e., the state residuals), and 
measurement errors. LST models have been previously applied to 
assess transsituaional consistency of psychological measures 
(Steyer, Schmitt, & Eid, 1999) and of biological measures, such 
as asymmetries in electroencephalography (Hagemann, Hewig, 
Seifert, Naumann, & Bartussek, 2005; Hagemann, Naumann, 
Thayer, & Bartussek, 2002), cerebral blood flow in resting state 
(Hermes et al., 2009), or salivary Cortisol awakening responses 
(Hellhammer et al., 2007; Kirschbaum et al., 1999). 

We measured electrocardiographs of 60 healthy student partici- 
pants at three measurement occasions once a week in three differ- 
ent measurement conditions: spontaneous breathing, paced 
breathing, and during a reaction time task. Thus, the measurement 
situation varied within each participant by occasion and experi- 
mental condition. On the basis of previous results (Pinna et al., 
2007; Sandercock et al., 2005), we hypothesized good reliabilities 
of HRV measures, that is, low portions of error variance. We also 
expected that besides trait and error variance, occasion- and 



condition-specific variance would significantly contribute to total 
HRV variance. 

Materials and Methods 

Participants 

Sixty healthy students (43 women, 17 men, M age = 22.2 years, 
SD = 1.9, range 19-27 years) of the University of Trier participated 
in the study, which consisted of three measurement occasions 
separated from each other by exactly 1 week. Participation was 
compensated with course credit. All participants underwent a lab- 
specific standardized screening interview prior to the experiment to 
check for exclusion criteria. In this interview, a trained psycholo- 
gist asked the participants (a) whether they currently suffered from 
any chronic or acute disease, (b) whether they currently took any 
medication other than oral contraceptives or occasional pain killers, 
and (c) whether they had ever suffered from any cardiovascular, 
metabolic, psychiatric, or stress-related complaints. The psycholo- 
gist also instructed all participants to refrain from smoking, 
caffeine, and meals 1 h prior to the experiment. According to rules 
of the local ethical committee, participants confirmed by signing 
that they had adequately and completely informed the experimenter 
of their medical status and that they will have to inform him 
immediately if their medical status may change. 

Procedure 

The experimental procedure was kept constant for each participant 
and each measurement occasion (Y1-Y3). After arriving in the 
laboratory, participants were prepared for the ECG measurement. 
Besides ECG electrodes, devices for electrocardiac impedance, 
breathing frequency, and blood pressure measurements were 
applied. The experiment consisted of three 6-min measurement 
conditions: (1) spontaneous breathing (SB), (2) paced breathing 
(PB), and (3) a reaction time task (RTT). Participants were told to 
sit still and relax during the spontaneous breathing condition. For 
the paced breathing condition, participants were instructed to 
breath in a 0.25-Hz rhythm via headphones (inspiration-expiration 
length ratio: 2:3). We included this condition because it is well 
known that the large within-subject variations in respiratory rates 
seen in spontaneous breathing conditions may confound the 
relationship between high frequency HRV measures and cardiac 
vagal tone (e.g., Grossman, Stemmler, & Meinhardt, 1990; Ritz & 
Dahme, 2006; Ritz, Thons, & Dahme, 2001). The reaction time 
task consisted of an adaptive choice reaction time task, which has 
been described earlier (Schachinger, Cox, Linder, Brody, & Keller, 
2003). Participants responded to the presentation of colored lights 
that appeared in random sequence by pressing corresponding 
buttons as accurately and quickly as possible. Using a PC-based 
control algorithm, the interstimulus intervals were shortened or 
lengthened, thereby modifying task difficulty so that subjects' false 
response rate within a continuously moving window approached 
50%. All experiments started with spontaneous breathing, followed 
by paced breathing, and ending with the reaction time task. 

Data Recording and Quantification 

Physiological data were recorded by Physiocorder software (Insti- 
tute of Biomedical Engineering, University of Stuttgart, Stuttgart, 
Germany) at a 1-kHz sampling rate and 16 bit resolution via 
a customized amplifier device (Curio Medizinelektronik, Bonn, 
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Germany) and stored to hard disk. ECG hardware band-pass filter 
settings were 0.5 to 500 Hz. Standard Ag/AgCl electrodes (ECG 
Tyco Healthcare H34SG Ag/AgCl electrodes of 45 mm dia- 
meter) were used for ECG measurements. ECG electrodes were 
attached according to a standard lead II configuration. Respiratory 
activity was measured via a thoracic respiratory belt (James Long 
Company, Inc., Caroga Lake, NY) placed between the fifth and 
eighth ribs. The signal was band-pass filtered (0.05-10 Hz) before 
it was digitized. Psychomotor responses during the reaction time 
task were recorded via a Serial Response Box device (Psychology 
Software Tools, Inc., Sharpsburg, PA), and responses were 
recorded by E-Prime software (see above). 

ECG data were processed with WinCPRS software (WinCPRS 
1.16, Absolute Aliens Oy, Turku, Finland). Interbeat intervals were 
calculated from the ECG and manually corrected, with a normal 
cycle RR-interval time series as output signal. Time domain meas- 
ures were directly calculated from the RR-interval series. The 
spectral analysis of the RR-interval series was done with WinCPRS 
software using a Fast Fourier Transform (FFT) routine. The 
RR-interval time series was linearly interpolated and resampled 
with a sampling rate of 5 Hz, the resampled data were tapered using 
a Hanning window, and the windowed data zero padded to the next 
power of 2. The FFT spectrum was smoothed using a sliding 
triangular weighting function in order to increase the number of 
degrees of freedom and thus improve the statistical relevance of 
the spectrum. 

We defined the high frequency band (HF) as 0.14 to 0.4 Hz 
and used the following time and frequency HRV parameters for 
statistical analyses: root mean square of successive differences 
(RMSSD), percent of difference between adjacent RR intervals that 
are greater than 50 ms (pNN50), and power in the high frequency 
band (HF power). We focused on those parameters because they are 
most commonly used to estimate parasympathetic control over 
heart rate (cf. Task Force, 1996). We used data of the whole 6-min 
segments of each measurement condition (SB, PB, RTT) at each 
occasion (Y1-Y3) for the calculation of ICC. For the structural 
equation modeling (SEM), we divided each 6-min segment into 
two 3-min segments as the estimation of latent state-trait models 
requires two identical measures at each measurement occasion 
(Steyer et al., 1999). This results in six measures for each of the 
three conditions (SB, PB, RTT) with Yn and Y 2] serving as two 
parallel measures in the first measurement occasion, Y ]2 and Y 22 as 
the two parallel measures in the second measurement occasion, 
and Yn and Y 23 as two parallel measures in the third measurement 
occasion. 

The respiratory signal was processed with WinCPRS software, 
as well. Inspiration and expiration phases were automatically 
detected by the WinCPRS algorithm and manually confirmed. We 
calculated mean respiratory frequency for each participant as final 
output measure. 

Statistical Analyses 

Respiratory activity. We calculated a 3x3 repeated- 
measurement analysis of variance (ANOVA) with the measurement 
occasion and measurement condition and the dependent variable 
mean respiratory activity to test for differences in the breathing 
frequency, which may potentially influence HRV parameters 
reflecting parasympathetic activity. Where appropriate, the 
Greenhouse-Geisser procedure was applied to correct for potential 
violations of the sphericity assumption. Effect sizes of significant 
results are reported as proportion of explained variance (r) 2 ). 



Intraclass correlation coefficients. To compare our results to 
those of previous studies (e.g., Pinna et al., 2007; Schroeder et al., 
2004; Sloan et al., 1995; for a review, see Sandercock et al., 2005) 
we computed the ICC for each HRV parameter and measurement 
condition as 

ICC = (BMS - WMS)/(BMS + WMS), 

where BMS is the between- subject mean square and WMS is the 
within-subject mean square. 

SEM analyses. We computed LST parameters to quantify latent 
components. In LST models, it is presumed that the observed 
variables are affected by the person, by the situation, and the 
person-situation interaction and by measurement errors. In the 
present study, different situations emerge on different measure- 
ment occasions and different experimental conditions. 

Model 1. The first model addressed the questions of the pro- 
portions of variance in HRV explained by measurement errors, a 
transsituationally stable trait, and occasion-specific state effects 
separately for the three measurement conditions. Therefore, the 
observed variables (Y, t ), that is, the HRV measures in the first and 
second 3 min (t) of the three measurement occasions (k), are 
decomposed in two stages: In a first step, the observed variables are 
decomposed into a latent state (S*) and into measurement errors 
(eft). In a second step, this latent state (St) is decomposed into a 
latent trait (T) and into occasion-specific state residuals (SRa). The 
former (T) represents individual differences that are consistent 
across the three measurement situations (in this case, measure- 
ment occasions), and the latter (SR t ) represent the effects of the 
situation and person-situation interaction on the latent state. This 
first model is displayed in Figure 1 . 




Figure 1. Model 1 : Latent state-trait model for two indicators (;') and 
three occasions of measurement (k). Y ;1 are manifest variables (HRV 
measures). e, t are latent errors. S ( are latent states. T is the latent trait, and 
SRi are latent state residuals. For simplicity, all factor loadings Bg and y, t 
are set equal to 1. 
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We tested this LST model separately for the three measurement 
conditions (SB, PB, and RTT) to confirm the hypothesized com- 
position of the HRV measures as a compound of a transsituation- 
ally consistent trait and occasion-specific components. 

We also performed two additional LST models to analyze 
influences of the three measurement conditions (SB, PB, and RTT) 
on HRV. 

Model 2. The second model investigated the effects of the three 
different conditions, SB, PB, and RTT, on HRV measures at one 
measurement occasion. The first measurement occasion was 
chosen as participants were naive to the measurement setup. This 
corresponds best to the many studies in which single short-term 
HRV measures are used as a biological trait for psychological or 
physiological functioning. Similar to the first model (see Figure 1), 
this second model decomposes the observed variables in two 
stages: In a first step, the observed variables are decomposed into 
a latent state (Sj) and into measurement errors (e,j). In a second 
step, this latent state (Sj) is decomposed into latent trait (T) and into 
condition-specific residuum (SR;). Only this time the observed 
variables (Y,,) represent HRV measures in the first and second 
3 min (;') in the three conditions (/') at the first occasion. Again, the 
trait (T) represents individual differences that are consistent 
across the three measurement conditions (SB, PB, RTT), and the 
condition-specific state residuals (SR ; ) represent the effects of the 



situation (in this case, measurement condition) and person- 
situation interaction on the latent state. The purpose of this model 
was to estimate the common trait variance across the three meas- 
urement conditions as well as the portions of variance specific to 
each of the measurement conditions. 

Model 3. In the third model, we analyzed the generalizability 
of stable individual differences across the three measurement con- 
ditions and measurement occasions. This model decomposed the 
manifest variables (Y,a), that is, HRV measures in the three condi- 
tions (SB, PB, RTT) at the three occasions, into a condition-specific 
latent trait (T,), occasion-specific state residuals (SR*), and meas- 
urement errors (e#). As shown in Figure 2, we allowed correlations 
between the three condition-specific trait variables (T/). These cor- 
relations indicate the generalizability of consistent individual dif- 
ferences across conditions (greater correlations suggest greater 
generalizability). The portion of occasion-specific state residuals 
indicates the generalizability of consistent individual differences 
across occasions (smaller variances of the occasion-specific state 
residuals suggest greater generalizability). Please note that the 
whole 6-min segments of each condition (SB, PB, RTT) at each 
measurement occasion were used to estimate this third model. 

The SEM analyses were based on the covariance matrices of the 
parallel HRV measures. We fitted the LST models to the covariance 
matrix by minimizing the generalized least squares (GLS) discrep- 
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Table 1. Means (M) ± one standard derivation (SD) as well as intraclass correlation coefficients (ICC) 
for heart rate variability parameters 



95% confidence interval ( CI) 



Measurement 1 (Ml ± SD\) Measurement 2 (M2 ± SD2) Measurement 3 (M3 ± SD3) 



ICC (95% CI) 



Spont. breath. 



Paced breath. 



Rt task 



ln-RMSSD 
ln-pNN50 
ln-HF power 
ln-RMSSD 
ln-pNN50 
ln-HF power 
ln-RMSSD 
ln-pNN50 
ln-HF power 



3.7 

3.3 : 

6.4 : 
3.6 : 
3.3 : 

6.2 : 

3.3 : 
2.9 : 
5.2 : 



0.4 
0.5 
0.9 
0.4 
0.5 
0.8 
0.5 
0.5 
1.0 



3.5 : 
3.2 : 

6.1 : 

3.6 : 

3.2 : 

6.1 : 

3.2 : 
2.8 : 
5.0 : 



0.5 
0.6 
1.0 
0.5 
0.5 
1.0 
0.5 
0.5 
1.1 



3.4 : 

3.0 : 
5.8 : 
3.4 : 

3.1 : 
5.8 : 
3.1 : 



0.5 
0.5 
0.9 
0.5 
0.5 
0.8 
0.8 



2.7 ± 0.5 

4.8 ± 1.0 



0.71 (0.54- 
0.73 (0.59- 
0.70 (0.53- 
0.72 (0.57- 
0.71 (0.56- 
0.73 (0.58- 
0.85 (0.75- 
0.83 (0.74- 
0.82 (0.72- 



0.82) 
0.83) 
0.81) 
0.82) 
0.82) 
0.83) 
0.90) 
0.90) 
0.89) 



Note. Spont. breath. = spontaneous breathing; paced breath. = paced breathing; rt task = reaction time task; RMSSD = root mean square of successive 
differences (ms); pNN50 = percent of difference between adjacent NN intervals that are greater than 50 ms (%); HF power = high frequency power (ms 2 ). 



ancy function. 1 Each model was analyzed in different versions, 
beginning with a maximally restricted model version (equal vari- 
ances of the measurement errors, equal variances of the state 
residuals, equal effects of states, equal effects of traits). If 
the model could not be accepted, a stepwise liberalization of these 
restrictions was introduced (Hagemann et al., 2002; Schmitt & 
Steyer, 1993). Note that these liberalizations are portions of the 
prespecified general model. No "a posteriori" model liberalizations 
were introduced to prevent nonreplicable and errorneous model 
solutions (MacCallum, 1995). We evaluated the general model fit 
with chi-square statistics and goodness-of-fit indices (comparative 
fit index [CFI], Bentler, 1990; root-mean- square error of approxi- 
mation [RMSEA], and the probability of the observed RMSEA 
under the null hypothesis RMSEA < .05 [pRMSEA]; Browne & 
Cudeck, 1993). The significance (a = .05) of single model param- 
eters was evaluated with the critical ratio (C.R.) of a parameter and 
its standard error (C.R. = Var/SE). 

On the basis of the accepted models, LST parameters were 
computed (for formulas, see Steyer et al., 1992; Steyer & Schmitt, 
1990). A coefficient of reliability specifies the proportion of vari- 
ance of the manifest variable that is due to all error-free latent 
components. A coefficient of trait specificity 2 represents the per- 
centage of variance of the manifest variable that is determined by 
transsituationally consistent individual differences. In contrast, a 
coefficient of occasion specificity reflects the percentage of vari- 
ance of the manifest variables that is determined by effects of the 
situation and the person-situation interaction in the first and third 
models, where the situations varied by occasion. In a similar vein, 
a coefficient of condition specificity reflects the percentage of 
variance of the manifest variables that is determined by effects of 
the situation and the person-situation interaction in the second 
model, where the situation varied by condition. Please note that we 
corrected the split-half in Models 1 and 2 analogously to the 
Spearman-Brown prophecy formula (see Steyer & Schmitt, 1990). 



1. We used the GLS rather than the maximum likelihood method 
because a simulation study showed that the GLS approach yields more 
precise results if the sample size is rather small and the multivariate nor- 
mality assumption does not hold (Hu & Bentler, 2005). Note that there is no 
strict test for the multivariate normality assumption and that the present 
sample is rather small. 

2. This coefficient is usually termed "common consistency" (e.g., 
Steyer et al., 1992). However, we followed a suggestion of Hagemann et al. 
(2002) and termed this coefficient "trait specificity" because it quantifies 
the proportion of the variance that is due to individual differences in a latent 
trait, whereas this trait is the factor that is transsituationally consistent and 
temporally stable across the measurement occasions and conditions. 



All statistical analyses were conducted with SPSS for Windows 
(Version 18, SPSS Inc.); SEM analyses were performed with Amos 
(Version 18, Amos Development Corporation). 

Results 

Mean Respiratory Activity 

Seven participants provided incomplete breathing data because of a 
technical malfunction or major artifacts and were thus excluded 
from analysis. Participants had a higher mean respiratory activity 
during the reaction time task (At = 0.36 Hz, SD = 0.08 Hz) 
F(2,104) = 97.71, p < .001, e = .92, if = .65, than during the two 
remaining conditions, which were (by chance) almost identical 
(spontaneous breathing: At = 0.25 Hz, SD = 0.06 Hz; paced breath- 
ing: At =0.25 Hz, 5D = 0.01Hz). All participants followed the 
instructions for inspiration and expiration in the paced breathing 
condition and thus had nearly identical breathing frequencies 
at all measurement occasions (Yi, Y 2 , and Y 3 : At = 0.25 Hz, 
SD = 0.01 Hz). A significant Occasion x Condition interaction, 
F(4,208) = 4.34, p = .018, £ = .47, if = .08, indicates a decreased 
breathing frequency during the first measurement occasion (Yi: 
M = 0.24 Hz, SD = 0.06 Hz vs. Y 2 : M = 0.25 Hz, SD = 0.06 Hz: 
p = .02) of the spontaneous breathing condition (Y 3 : At = 0.25 Hz, 
SD = 0.06 Hz). Furthermore, the mean respiratory frequencies 
during the reaction time task condition (Y^ M = 0.37 Hz, SD = 
0. 1 1 Hz; Y 2 : At = 0.36 Hz, SD = 0.07 Hz; Y 3 : At = 0.34 Hz, SD = 
0.06 Hz) did not differ over the three measurement occasions 
(ps> .10). 

Descriptive Statistics and ICCs 

All HRV parameters had a markedly right-skewed distribution 
(Shapiro-Wilk's p < .05) and were thus log-transformed (In). 
Means and standard derivations as well as ICC coefficients of the 
log-transformed HRV parameters in the three conditions at the 
three measurement occasions are presented in Table 1. ICCs 
were between 0.70 and 0.85 and are thus comparable to those re- 
ported in previous studies (e.g., Pinna et al., 2007; Schroeder et al., 
2004; Sloan et al., 1995; for a review, see Sandercock et al., 2005). 

LST Models of HRV 

Model 1. The purpose of the first model (see Figure 1) was to 
estimate the portion of variance explained by individual differences 
in a transsituationally consistent trait and occasion-specific state as 
well as the reliability of HRV measures. 
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Table 2. Model 1: Chi-Square Statistics, Goodness-of-Fit Indices, and Estimated Model Parameters of a Latent State-Trait Model of 
Heart Rate Variability for Three Parameters in Three Measurement Conditions 



Goodness of fit Estimated model parameters 



Task 


Measure 


X 2 (18,/V=60) 


P 


CFI 


RMSEA 


Prmsea 


Var(y) 


Var(e) 


Var(T) 


Var(SR) 


Spont. breath. 


ln-RMSSD 


19.11 


.385 


0.99 


0.03 


.551 


0.164 


0.010 


0.085 


0.078 




SE 














0.001 


0.023 


0.012 




ln-pNN50 


10.49 


.9 1 5 


1.00 


0.00 


.958 


0.275 


0.017 


0.136 


0.122 




SE 














0.002 


0.034 


0.017 




ln-HF power 


10.10 


.928 


1.00 


0.00 


.965 


0.82 


0.100 


0.389 


0.328 




SE 














0.011 


0.101 


0.052 


Paced breath. 


ln-RMSSD 


26.07 


.098 


0.90 


0.09 


.202 


0.164 


0.011 


0.073 


0.080 




SE 














0.001 


0.021 


0.012 




ln-pNN50 


24.39 


.143 


0.92 


0.08 


.269 


0.225 


0.017 


0.095 


0.113 




SE 














0.002 


0.029 


0.017 




ln-HF power 


18.16 


.445 


1.00 


0.01 


.609 


0.654 


0.101 


0.339 


0.214 




SE 














0.011 


0.082 


0.039 


Rt task 


ln-RMSSD" 


24.36 


.110 


0.87 


0.09 


.217 


0.211/0.211/0.273 


0.020 


0.155 


0.036/0.036/0.098 




SE 














0.002 


0.034 


0.011/0.011/0.024 




ln-pNN50 b 


23.88 


.067 


0.84 


0.10 


.143 


0.177/0.166/0.213 


0.024/0.013/0.026 


0.120 


0.033/0.033/0.067 




SE 














0.005/0.003/0.005 


0.029 


0.010/0.010/0.020 




ln-HF power 


16.71 


.54 


1.00 


0.00 


.70 


0.95 


0.18 


0.59 


0.18 




SE 














0.182 


0.134 


0.181 



Note. All reported statistics are based on parametric generalized least squares estimates. The same restrictive latent state-trait model was analyzed for the HRV parameters of each 
condition (equal variances of the measurement error e,*, equal variances of the state residuals SRt, equal effects of states Si, equal effects yj, of trait T) with df- 18 and /V— 60, 
except for *df = 17, N — 60, one state residuum (SRa) variance was set free; b df= 15 N — 60, one state residuum (SRa) and one error variance (e,i) were specific of each measurement 
occasion. Therefore, separate coefficients of true score variance Var(y,A) 1,2 , state residuum Var(SRi) 1,2 (and error variance Var(e,i) 2 ) for the three measurement occasions have to 
be reported for these models. CFI = comparative fit index; RMSEA — root-mean-square error of approximation; /?rmsea = probability of the observed RMSEA under the null 
hypothesis RMSEA < .05. 



The restrictive LST model (equal variances of the measurement 
errors and state residuals, equal effects of states and of trait) showed 
acceptable fit to the data for all parameters in all measurement 
conditions, all % 2 s(18, N = 60) < 26.07, ps > .10, CFIs > .90, 
RMSEAs < .09, prmseaS > .20) except for RMSSD and pNN50 in 
the reaction time task. An LST model with one liberalized state 
residuum (SRt) showed acceptable fit for RMSSD in the reaction 
time task condition, x 2 (ll, N = 60) = 24.36, p=.U, CFI = .87, 
RMSEA = .09, Prmsea = -22. A liberalization of one state residuum 
(SRt) and the acceptance of measurement occasion-specific error 
variance (e») were necessary to get an acceptable fit for pNN50 in 
the same condition, x 2 (15, N = 60) = 23.88, p = .07, CFI = .84, 
RMSEA = .10, Prmsea = -14. Table 2 summarizes the respective 
model test statistics for all parameters and conditions and addition- 
ally presents estimated variances and their standard errors for the 
latent measurement errors, state residuals, and traits. For all param- 
eters and conditions, the latent error variances, the latent state 



residual variances, and the latent trait variances contributed signifi- 
cantly to the variance of the manifest measures (all C.R.s & 3.29). 

Reliability, occasion specificity, and trait specificity. The 
parameter estimates for the LST model (see Table 2) allowed the 
computation of coefficients of reliability, occasion specificity, and 
trait specificity, which are presented in Table 3. Reliability was in 
the range of .81 to .99 (Mdn = .92) for all HRV parameters and 
measurement conditions. The occasion specificity was in the range 
of .17 to .50 (Mdn = 0.40), and the trait specificity was in the range 
of .42 to .73 (Mdn = 0.52) for all parameters. Interestingly, esti- 
mated occasion specificity was lowest and trait specificity was 
largest in the reaction time task for all HRV parameters (occasion 
specificity range .17-36, Mdn = 0.23; trait specificity range .56- 
.73, Mdn = .65). Thus, up to 50% of individual differences in a 
single measurement of HRV could be explained by situation and 
person-situation interaction. 



Table 3. Model 1: Reliability, Expected Occasion Specificity and Trait Specificity of HRV Measures after an Aggregation of the Data 
across Two and Three Measurement Occasions 



One occasion Two occasions Three occasions 



Task 


Measure 


Rel(y) 


OccSpec(y) 


TraSpe(y) 


OccSpec(y) 


TraSpe(y) 


OccSpec(y) 


TraSpe(y) 


Spont. breath. 


ln-RMSSD 


0.994 


0.476 


0.518 


0.314 


0.659 


0.234 


0.743 




ln-pNN50 


0.938 


0.444 


0.495 


0.306 


0.662 


0.233 


0.746 




ln-HF power 


0.878 


0.401 


0.476 


0.290 


0.645 


0.224 


0.732 


Paced breath. 


ln-RMSSD 


0.933 


0.488 


0.445 


0.349 


0.616 


0.270 


0.706 




ln-pNN50 


0.924 


0.502 


0.422 


0.367 


0.594 


0.287 


0.687 




ln-HF power 


0.846 


0.327 


0.518 


0.234 


0.683 


0.179 


0.764 


Rt task 


ln-RMSSD" 


0.905/0.905/0.927 


0.171/0.171/0.359 


0.735/0.735/0.568 


0.103/0.103/0.238 


0.847/0.847/0.724 


0.074/0.074/0.177 


0.893/0.893/0.798 




ln-pNN50 b 


0.864/0.922/0.878 


0.186/0.199/0.315 


0.678/0.724/0.563 


0.119/0.120/0.214 


0.808/0.839/0.721 


0.087/0.086/0.161 


0.863/0.887/0.795 




ln-HF power 


0.809 


0.190 


0.619 


0.130 


0.764 


0.100 


0.830 



Note. The coefficients of occasion specificity, OccSpe(Y), and trait specificity, TraSpe(Y), of the aggregated HRV measures are based on the estimated model parameters of the 
restrictive latent state-trait model for HRV measures except for ;l and b and were computed following Steyer and Schmitt (1990). As in a one state residuum (SR*) variance and 
in b one state residuum (SRt) variance and one error variance (e&) are specific of each measurement occasion, coefficients of reliability, occasion and trait specificity have to be 
reported separately for each measurement occasion. N - 60. 
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Aggregation across occasions. The coefficients of occasion 
specificity and trait specificity were recomputed for the all HRV 
measures for the case that these variables are aggregated across two 
and three occasions of measurement. The resulting coefficients are 
presented in Table 3. After an aggregation across two occasions, 
the expected occasion specificity was in the range of .10 to .37 
(Mdn = .29) across all parameters; after an aggregation across three 
occasions it was in the range of .07 to .29 (Mdn = .22). After an 
aggregation across two and three occasions, however, the trait 
specificity was in the range of .59 to .85 (Mdn = .66) and .71 to .89 
(Mdn = .75), respectively. Hence, the portions of variance in HRV 
measures explained by transsituationally stable individual differ- 
ences (trait) increases from 52% at one single measurement to 
about 66% or 75% when aggregated over two or three measure- 
ment occasions. 

Model 2. This model estimated the portion of condition-specific 
variances and trait variance of HRV measures at a single — the 
first — measurement occasion. The restrictive LST model (equal 
variances of the measurement errors and state residuals, equal 
effects of states and of trait) showed acceptable fit to the data 
for HF power, % 2 (18, N = 60) = 17.87, p = .40, CFI = .98, 
RMSEA = .03, prmsea = .56, and pNN50, % 2 (18, N= 60) = 17.93, 
p = .33, CFI = .97, RMSEA = .05, prmsea = .49. An LST model 
with condition-specific measurement errors showed acceptable fit 
to the data for RMSSD, % 2 (14, N = 60) = 20.92, p = .10, CFI = .88, 
RMSEA = .09, Prmsea = -20. Table 4 presents estimated variances 



Table 4. Model 2: Estimated Model Parameters of a Latent 
State-Trait Model of Heart Rate Variability for Three Parameters 
in the Three Measurement Conditions ('$) Spontaneous Breathing 
(Spont. Breath., SB), Paced Breathing (Paced Breath., PB), and 
Reaction Time Task (Rt Task, RTT) at the First Measurement 
Occasion 



Estimated model parameters 



Task 


Measure 


Var(y) 


Var(T) 


Var(SR) 


Var(e) 


Spont. breath. 


ln-RMSSD" 


0.127 


0.093 


0.027 


0.007 




SE 




0.023 


0.010 


0.002 




ln-pNN50 


0.196 


0.120 


0.056 


0.020 




SE 




0.031 


0.019 


0.002 




ln-HF power 


0.778 


0.448 


0.209 


0.121 




SE 




0.094 


0.063 


0.014 


Paced breath. 


ln-RMSSD" 


0.124 


0.093 


0.018 


0.013 




SE 




0.023 


0.009 


0.003 




ln-pNN50 


0.173 


0.120 


0.033 


0.020 




SE 




0.031 


0.016 


0.002 




ln-HF power 


0.569 


0.448 


0.000 


0.121 




SE 




0.094 


0.000 


0.014 


Rt task 


ln-RMSSD a 


0.175 


0.093 


0.054 


0.028 




SE 




0.023 


0.017 


0.006 




ln-pNN50 


0.220 


0.120 


0.080 


0.020 




SE 




0.031 


0.022 


0.002 




ln-HF power 


0.998 


0.448 


0.429 


0.121 




SE 




0.094 


0.103 


0.014 



Note. All reported statistics are based on parametric generalized least 
squares estimates. The same restrictive latent state-trait model was ana- 
lyzed for the HRV parameters (equal variances of the measurement error e,j, 
equal variances of the state residuals SR,, equal effects B,y of states S, equal 
effects yj of trait T) with df= 18 and N = 60, except for a rf/ = 14, N= 60 
where the error variance (e,j) was specific for the measurement condition. 
Therefore, specific coefficients of error variance 11 for the three measure- 
ment conditions are reported for this model. 



Table 5. Model 2: Reliability, Condition Specificity, and Trait 
Specificity of HRV Measures at the First Measurement Occasion 



Toe I' 


Measure 


Kelty) 


ConSpe(y) 


IraSpe(y) 


Spont. breath. 






0 911 


U. / jL 




ln-pNN50 


0.898 


0.286 


0.612 




ln-HF power 


0.844 


0.269 


0.576 


Paced breath. 


ln-RMSSD 


0.895 


0.145 


0.750 




ln-pNN50 


0.884 


0.191 


0.694 




ln-HF power 


0.787 


0.000 


0.787 


Rt task 


ln-RMSSD 


0.840 


0.309 


0.531 




ln-pNN50 


0.909 


0.364 


0.545 




ln-HF power 


0.879 


0.430 


0.449 



Note. The coefficients of condition specificity, ConSpe(Y), trait specificity, 
TraSpe(Y), of the HRV measures are based on the estimated model param- 
eters of the latent state-trait model for HRV measures described in Table 4 
and were computed following Steyer and Schmitt (1990). N = 60. 



and their standard errors for the latent measurement errors, state 
residuals, and traits. For all parameters, the latent error variances, 
the latent state residual variances, and the latent trait variances 
contributed significantly to the variance of the manifest measures 
(all C.R.s > 2.00), except for the state residual of the paced breath- 
ing condition in HF power (SRpb, C.R. = 1.25, p = .21). Therefore, 
the model was reestimated after setting the variance of SRpb to 0. 

Reliability, condition specificity, and trait specificity. We also 
computed coefficients of reliability, condition specificity, and trait 
specificity with the parameters estimated for this second LST 
model (see Table 5). Reliability was in the range of .79 to .95 
(Mdn = .90) for all HRV parameters across the three measurement 
conditions. Condition specificity was in the range of 000 to .36 
(Mdn = .27) and thus explained about 27% of the interindividual 
variance in HRV. Condition-specific effects were greater in the 
reaction time task (Mdn = .40) than in spontaneous (Mdn = .20) 
and paced breathing (Mdn = .23). Trait specificity (i.e., the variance 
proportion that is due to a factor that is consistent across the three 
measurement conditions) was in the range of .45 to .79 (Mdn = .61) 
across all HRV parameters and conditions and thus explained about 
61% of individual differences in HRV. Trait effects were greater for 
spontaneous (Mdn = .70) and paced (Mdn = .69) breathing than in 
the reaction time task (Mdn = .55). Thus, about 61% of interindi- 
vidual differences in a single measurement of HRV could be 
explained by condition-unspecific trait variance, whereas about 
27% of variance was explained by the measurement condition and 
person-condition interactions. 

Model 3. The third model estimated the generalizability of stable 
individual differences across the three measurement conditions 
and measurement occasions. An LST model with equal variances 
of the measurement errors, condition-specific trait variances, 
and occasion-specific state variances showed acceptable fit for 
all parameters: ln-RMSSD, x 2 (35, N= 60) = 34.70, p=.48, 
CFI =1.00, RMSEA = .00, prmsea = .70; ln-pNN50, x 2 (35, 
N= 60) = 43.08, p = .16, CFI = 90, RMSEA = .06, /?rmsea=.35; 
ln-HF power, % 2 (35, N = 60) = 33.90, p = .52, CFI =1.00, 
RMSEA = .00, Prmsea = -73. Table 6 presents estimated variances 
and their standard errors for the latent measurement errors, state 
residuals, and traits. For all parameters, the latent error variances, 
the latent state residual variances, and the latent trait variances 
contributed significantly to the variance of the manifest measures 
(all CRs > 2.48). 
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Table 6. Model 3: Estimated Parameters of a Latent State-Trait Model of Heart Rate Variability for Three Parameters in the Three 
Measurement Conditions Spontaneous Breathing (Spont. Breath., SB), Paced Breathing (Paced Breath., PB), and Reaction Time Task (Rt 
Task, RTT) at the Three Measurement Occasions Yl, Y2, and Y3 



Estimated model parameters 



Occasion 


Measure 


Var(y SB ) 


Var(y PB ) 


Var(y RTT ) 


Var(T SB ) 


Var(TpB) 


Var(T RT T) 


Var(SB) 


Var(e) 


Yl 


ln-RMSSD 


0.146 


0.144 


0.209 


0.093 


0.091 


0.156 


0.040 


0.013 




SE 








0.023 


0.022 


0.035 


0.015 


0.001 




In-pJNJN jU 


0.222 


0.224 


0.215 


0.135 


0.137 


0.128 


0.06/ 


0.020 




SE 








0.033 


0.03 1 


0.030 


0.020 


0.002 




ln-HF power 


0.691 


0.577 


0.912 


0.437 


0.323 


0.658 


0.138 


0. 1 16 




SE 








0.101 


0.083 


0.142 


0.056 


0.013 


Y2 


ln-RMSSD 


0.181 


0.179 


0.244 


0.093 


0.091 


0.156 


0.075 


0.013 




SE 








0.023 


0.022 


0.035 


0.020 


0.001 




ln-pNN50 


0.215 


0.217 


0.208 


0.135 


0.137 


0.128 


0.060 


0.020 




SE 








0.033 


0.031 


0.030 


0.020 


0.002 




ln-HF power 


0.759 


0.645 


0.980 


0.437 


0.323 


0.658 


0.206 


0.116 




SE 








0.101 


0.083 


0.142 


0.071 


0.013 


Y3 


ln-RMSSD 


0.187 


0.185 


0.250 


0.093 


0.091 


0.156 


0.081 


0.013 




SE 








0.023 


0.022 


0.035 


0.021 


0.001 




ln-pNN50 


0.201 


0.203 


0.194 


0.135 


0.137 


0.128 


0.046 


0.020 




SE 








0.033 


0.031 


0.030 


0.020 


0.002 




ln-HF power 


0.753 


0.639 


0.974 


0.437 


0.323 


0.658 


0.200 


0.116 




SE 








0.101 


0.083 


0.142 


0.064 


0.013 



Reliability, occasion specificity, trait specificity, and correla- 
tions between traits. Table 7 summarizes coefficients of reliability, 
occasion, and trait specificity computed with the parameters esti- 
mated for the third LST model. Reliability was in the range of .79 
to .95 (Mdn = .91), occasion specificity was in the range of .19 to 
.43 (Mdn = .29), and trait specificity was between .49 and .74 
(Mdn = .62). Thus, between 19% and 43% of interindividual dif- 
ferences in the three HRV parameters was explained by interindi- 
vidual differences in effects of the measurement occasion or 
person-occasion interactions. In addition, between 49% and 74% 
of individual differences in HRV parameters were explained by 
interindividual differences in condition-specific trait variance 
across the three measurement occasions. The three condition- 
specific traits were significantly intercorrelated (.52 < r < .92) for 
all HRV parameters. All parameters showed strongest correlations 
between spontaneous and paced breathing (.77 < r < .92). This 
indicates that the generalizability of stable individual differences 
across conditions is high — in particular regarding spontaneous and 
paced breathing conditions. 

Discussion 

The present study demonstrated good to excellent reliability (.81- 
.99) for parasympathetic HRV measures over the course of 3 weeks 



in healthy student participants. Besides very low estimated error 
variance, LST models also revealed that about one third to one half 
of the variance of a single HRV measurement can be explained by 
occasion-specific effects. This significant portion of occasion- 
specific variance may at least partly explain the mixed results with 
regard to the temporal stability of HRV measures as well as to the 
relationship between HRV measures and physiological and psycho- 
logical functions. Repeated measurements in well-controlled 
situations may be a possible recommendation for researchers inter- 
ested in these variables. 

One aim of the study was to estimate the temporal stability of 
HRV measures. We calculated intraclass correlation coefficients 
and found sufficiently great values for all investigated HRV 
parameters — RMSSD, pNN50, and HF power — over the course of 
3 weeks in three different task conditions. The ICCs of the present 
study showed better stability for all HRV parameters in the reaction 
time task (.82-85) than in the spontaneous (.70-.73) or paced 
breathing conditions (.7 1— .73). Although our results on spontane- 
ous and paced breathing were comparable to findings of earlier 
studies (e.g., Pitzalis et al., 1996; Schroeder et al., 2004; Sloan 
et al., 1995), the increased stability of HRV during task perform- 
ance was contrary to findings of two earlier studies (Li et al., 2009; 
Sloan et al., 1995). Both reported lower stability of high frequency 
HRV (9 months; 1.5 years) in response to different mental stress 



Table 7. Model 3: Reliability, Occasion Specificity, and Trait Specificity of HRV Measures 



Occasion 


Measure 


Rel(ysu) 


Rel(y PB ) 


Rel(y RTT ) 


OccSpe(y SB ) 


OccSpe(y PB ) 


OccSpe(y RTT ) 


TraSpe(y SB ) 


TraSpe(y PB ) 


TraSpe(y RTT ) 


Yl 


ln-RMSSD 


0.911 


0.910 


0.938 


0.274 


0.278 


0.191 


0.637 


0.632 


0.746 




ln-pNN5() 


0.910 


0.911 


0.907 


0.302 


0.299 


0.312 


0.608 


0.612 


0.595 




ln-HF power 


0.832 


0.799 


0.873 


0.200 


0.239 


0.151 


0.632 


0.560 


0.721 


Y2 


ln-RMSSD 


0.928 


0.927 


0.947 


0.414 


0.419 


0.307 


0.514 


0.508 


0.639 




ln-pNN5() 


0.907 


0.908 


0.904 


0.279 


0.276 


0.288 


0.628 


0.631 


0.615 




ln-HF power 


0.847 


0.820 


0.882 


0.271 


0.319 


0.210 


0.576 


0.501 


0.671 


Y3 


ln-RMSSD 


0.930 


0.930 


0.948 


0.433 


0.438 


0.324 


0.497 


0.492 


0.624 




ln-pNN5() 


0.900 


0.901 


0.897 


0.229 


0.227 


0.237 


0.672 


0.675 


0.660 




ln-HF power 


0.846 


0.818 


0.881 


0.266 


0.313 


0.205 


0.580 


0.505 


0.676 



Note. The coefficients of occasion specificity, OccSpe(Y), and trait specificity, TraSpe(Y), of the HRV measures are based on the estimated model parameters 
of the latent state-trait model of HRV measures described in Table 6 and were computed following Steyer and Schmitt (1990). N= 60. 
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tasks (reaction time task; arithmetic task; video game) than during 
rest. A possible explanation for these diverging results in the 
present study is that the "choice reaction time task" was a self- 
adaptive and well-controlled, cognitively demanding task, which 
may have guaranteed a similar performance and sympathetic/ 
parasympathetic activation pattern across the three measurement 
occasions. 

These conclusions are also consistent with the results from the 
LST analyses, which revealed good reliability for all task condi- 
tions, but larger proportions of estimated occasion specificity for 
spontaneous and paced breathing than for the reaction time task 
condition. Thus, occasion-specific influences (i.e., situation and 
person-situation interaction) had a stronger impact on HRV meas- 
ures in the less controlled task conditions. This may explain the 
substantial differences in reliability estimations via ICC found in 
the literature (for review, see Sandercock et al., 2005). Indeed, in 
the classical measurement model, test-retest temporal variability is 
included in the error term, and the ICC reflects only its relative 
magnitude. Pinna et al. (2007), for instance, reported ICCs of .70 to 
.88 but that increases as great as about 3.5 times and decreases of 
as much as 30% may occur between two HRV measurements taken 
1 day apart from each other. 

To our knowledge, the current study is the first to estimate the 
proportions of situation-specific factors and of consistent physi- 
ological traits using latent-state-trait equation modeling. In par- 
ticular, about 40% of the variance of HRV measures was due to 
occasion-specific effects, which suggests a considerable impact of 
the measurement situation and person-situation interactions on 
total HRV variance. The aggregation over two and three measure- 
ment occasions increased the portion of transsituationally consist- 
ent trait variance from about 49% for a single measurement to 66% 
and 75% across two and three occasions, respectively. The present 
data may therefore raise questions about the psychometric proper- 
ties of a single HRV measurement for assessing an individual's 
transsituationally consistent adaptability to situations. Our results 
suggest that a single measurement of an RR time series reflects a 
large amount of state variance. Therefore, we suggest the assess- 
ment of and aggregation over at least two measurements for the 
interpretation of high frequency HRV measures as a transsituation- 
ally stable biological marker or trait variable. Similar recommen- 
dations already exist for the assessment of tonic Cortisol levels 
(Cortisol awakening responses) with salivary Cortisol samples 
(Hellhammer et al., 2007) or for lateralization in tonic cortical 
activity with electroencephalography (Hagemann et al., 2002, 
2005). 

We also performed additional analyses to estimate influences of 
condition-specific variances on HRV. These models revealed that 
besides trait and error variance, occasion-specific variance contrib- 
uted significantly to total HRV variance. The second model showed 
that condition-specific variances accounted for about 27% of inter- 
individual differences in parasympathetic HRV measures whereas 
individual differences in condition-unspecific trait variance 
explained about 61% of interindividual differences. Again, the 
portion of variance explained by measurement error was small, 
and, in turn, the estimated reliability was good. Thus, even across 
the three difference measurement conditions, spontaneous breath- 
ing, paced breathing, and reaction time task, latent condition- 
unspecific trait variance explained almost two thirds of the variance 
of HRV. 

Finally a third model estimated the portion of interindividual 
differences in parasympathetic HRV measures explained by 
(occasion-specific) state variance and measurement condition- 



specific consistent trait variance. The results of this model revealed 
for all of the three HRV parameters that both interindividual dif- 
ferences in the effects of the particular situation at each measure- 
ment occasion as well as interindividual differences in stable traits 
had significant influences on the interindividual differences in HRV 
measures. The latter were specific for each measurement condition 
(spontaneous breathing, paced breathing, and reaction time task) 
and explained between 49% and 74% of the interindividual differ- 
ences in HRV measures. The state variance represented effects 
specific to the particular measurement occasion and explained 
between 19% and 43% of the interindividual differences in HRV 
measures. Most interesting were the high correlations between the 
three condition-specific traits, which indicate a generalizability of 
stable trait effects across the three measurement conditions. 

When interpreting these findings, one should bear in mind that 
in each of these alternative models the latent traits are defined as the 
common factors of different sets of manifest variables, which 
implies that the latent traits have different meanings. In particular, 
in the first model the latent trait is defined as the common factor of 
two parallel measures of the same construct across three occasions 
of measurement (the identical experimental condition was used in 
each measurement occasion). In the second model, the latent trait is 
defined as the common factor of two parallel measures across three 
experimental conditions (all taking place within the first measure- 
ment occasion). In the third model, each latent trait is defined as the 
common factor of one measure across three occasions of measure- 
ment (the identical experimental condition is used in each meas- 
urement occasion for each trait). Therefore, the meaning of the trait 
specificity and the situation (occasion or condition) specificity is 
not identical across the three models, which explains why the 
respective coefficients of specificity vary across the three models. 
Nonetheless, the total of the findings suggest that (a) about 30%- 
40% of variance of the HRV measures is due to effects of the 
situation and person-situation interaction (no matter whether the 
situation varies by measurement occasion or experimental condi- 
tion), (b) about 50%-60% of the variance of the HRV measures is 
due to individual differences in a latent trait, (c) the HRV measures 
show reliability estimates of about .90, and (d) the latent HRV traits 
show great correlations between experimental conditions, the latter 
of which may suggest that the latent HRV trait may be generalized 
across experimental conditions. 

Our results imply that situational and personal differences may 
substantially influence HRV measures and correspond well to those 
reported in behavioral genetic studies. These indicate that up to 
65% of the variance in HRV at rest and under stress can be 
explained by genetic influences (e.g., Kupper et al., 2004; Snieder, 
Boomsma, Van Doornen, & De Geus, 1997; Wang et al., 2009). 
The participants' current affective state could be an explanation for 
situational influences, as negative affective states may reduce high 
frequency HRV (Filaire, Protier, Massart, Ramat, & Teixeira, 2010; 
Schulz, Alpers, & Hofmann, 2008). Several researchers reported 
that participants may habituate to experimental settings (e.g., 
Langewitz, Riiddel, Noack, & Wachtarz, 1989), which could be 
interpreted as systematic influence. Furthermore, the menstrual 
cycle phase could be another systematic interindividual influence, 
which may have substantial impact on HRV (Bai, Li, Zhou, & Li, 
2009; Leicht, Hirning, & Allen, 2003; Sato, Miyake, Akatsu, & 
Kumashiro, 1995) but was not controlled in the present study. We 
suggest that this issue should be addressed in future research. In 
addition, other factors, such as respiratory rate, tidal volume, or 
momentary physical activity may influence a reliable estimation of 
cardiac vagal tone (for an extensive discussion, see Grossman & 
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Kollai, 1993; Grossman & Taylor, 2007). However, the here 
observed values of high reliability imply that these factors play a 
minor role in our data. 

The respiratory activities in the spontaneous and paced breathing 
conditions were almost identical, whereas it was increased in the 
reaction time task condition. Therefore, differences in absolute 
values or in temporal stability between the reaction time task and 
both resting conditions may also be influenced by the difference in 
the breathing frequency. Moreover, we did not observe any differ- 
ence in the respiratory activity between the three measurement 
occasions in each condition (0.25 Hz), except for lower mean res- 
piratory activity during the very first session (0.24 Hz). This differ- 
ence may, however, be associated with the above described 
habituation effect to laboratory experimental settings, which can 
often be observed in multiple laboratory measurements (e.g., Lange- 
witz et al., 1989). Finally, the similar average respiratory rates and 
average magnitudes of high frequency HRV measures in the spon- 
taneous and paced breathing conditions provide further evidence 
that paced breathing has little effect on respiratory sinus arrhythmia 
(RSA) and supports the use of pacing as a method to control for the 
sizable effects of breathing frequency on RSA that can severely 
confound estimation of cardiac vagal control of heart rate (see, e.g., 
Grossman et al., 1990; Ritz & Dahme, 2006; Ritz et al., 2001). 



Some limitations of the study need to be acknowledged. Effects 
of the experimental conditions should be interpreted with caution, 
as they were applied in a fixed order. This fixed order was neces- 
sary, because the validity of the spontaneous or paced breathing 
session as a "resting" condition would have been limited if they 
were performed after the reaction time task, which was considered 
as a stress-inducing task. However, revealing systematic differ- 
ences between the three conditions was not the aim of this study. 
Furthermore, parameters of temporal stability of parasympathetic 
HRV measures were based on three occasions with an interval of 1 
week. Therefore, further studies should address the issue of 
whether the temporal stability of HRV measures is still comparable 
if the interval between measurement occasions is larger. Finally, all 
drawn conclusions are limited to HRV parameters reflecting para- 
sympathetic activity. The temporal stability of other HRV measures 
was not addressed in this study. 

Taken together, these results suggest that HRV reflecting para- 
sympathetic activity is determined not only by trait-like factors and 
measurement error, but also by effects of the situation or the 
person-situation interaction. We therefore suggest an aggregation 
across at least two measurements when using parasympathetic 
HRV as a stable trait-like biomarker for cardiovascular health as 
well as emotional or cognitive functions. 
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